arXiv:1506.04994vl [astro-ph.SR] 16Jun2015 


Connecting Flares and Transient Mass Loss Events in 
Magnetically Active Stars 


Rachel A. Osten 1 
Space Telescope Science Institute 
3700 San Martin Drive, Baltimore, MD 21218 


ostenOstsci.edu 

and 

Scott J. Wolk 

Harvard-Smithsonian Center for Astrophysics 
60 Garden Street, Cambridge MA 02138 


Received 


accepted 


Also at Center for Astrophysical Sciences, Johns Hopkins University, Baltimore, MD 


21218 





ABSTRACT 


We explore the ramification of associating the energetics of extreme mag¬ 
netic reconnection events with transient mass loss in a stellar analogy with solar 
eruptive events. We establish energy partitions relative to the total bolomet- 
ric radiated flare energy for different observed components of stellar flares, and 
show that there is rough agreement for these values with solar flares. We apply an 
equipartition between the bolometric radiated flare energy and kinetic energy in 
an accompanying mass ejection, seen in solar eruptive events and expected from 
reconnection. This allows an integrated flare rate in a particular waveband to be 
used to estimate the amount of associated transient mass loss. This approach is 
supported by a good correspondence between observational flare signatures on 
high flaring rate stars and the Sun, which suggests a common physical origin. 
If the frequent and extreme flares that young solar-like stars and low-mass stars 
experience are accompanied by transient mass loss in the form of coronal mass 
ejections, then the cumulative effect of this mass loss could be large. We find 
that for young solar-like stars and active M dwarfs, the total mass lost due to 
transient magnetic eruptions could have significant impacts on disk evolution, 
and thus planet formation, and also exoplanet habitability. 


Subject headings: stars: activity — stars: flare — stars: late-type — stars: mass-loss 



3 


Introduction 


Changes in the magnetic configuration of the Sun’s outer atmosphere occur on short 
timescales and lead to the release of free energy which powers particle acceleration, plasma 
heating, shocks and mass motions. These solar eruptive events manifest themselves as 
flares, coronal mass ejections (CMEs), and proton events. Flares have been observed on 


many different ki nds of cool stars, ranging from F dwarfs (IMullan fe Mathioudakis 


G and K giants ([Ayres et ah 


flQsten et al. 


1999 


2001 


200411 . dM 


star-d isk interactions (Favata et al. 


limit (jStelzer et ah 


e flare stars (Haw 


Testa et ah 


ev et ah 


20001 ) . 


2007 ). tidally locked binary systems 


199511 . hyperactive young Suns containing 


200511 . and very low mass stars near the substellar 
2006h . Magnetic re connection processes are thought to unde rly the 


events, both on cool stars and the Sun ( Sh ib at a 


1999 


Shi bata fe Yoko vama 


1999 ). Detailed 


studies of various aspects of stellar flare phenomena show global agreement with processes 
occurring during solar flares, implying a similar physical origin, despite the disparity in 
scale - the largest stellar fl ares can be up to one million times more energetic than the 


largest solar flares (lOsten et ah 


white light stellar flares (Hawley et ah 


2007). Stell ar counterparts to solar flare phenom ena include 


region and_c 


iromospheric flares ( Hawley et ah 


(Osten et ah 

2005 

), nonthermal 

emission ( 

Osten & Bastian 

2006) 


1995). soft X-ray flares (iReale e t al. 19971), transition 


2003), nonthermal 


rard X-ray emission ( Osten et ah 


gyros ynchrotron emission 


20071), an d coherent radio 


2006 ). This diversity underscores the multi-wavelength nature 


of the flare process; the details connecting the energetics in the different wavelength regions 
is still unclear. 

While flares on cool stars are observed across the electromagnetic spectrum, the 
detection of transient mass loss signatur es in cool dwarf stars fa s with stea d y ma ss loss 
signatures) remains a difficult prospect (ILeitzinger et al.l 120111) . iHoudebinel (119961) and 


references therein list a number of outflows detected in stellar flares, but there are more 
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redshifted than blue-shifted emission components, and the hetergeneous nature of the 
events prevents a systematic exploration. Thus, in the absence of conclusive evidence of 
stellar CMEs, stellar astronomers must look to the Sun to gain insight into transient stellar 
mass loss on magnetically active stars. 


The commonality of stellar flare phenomenology with solar observations suggests 
that an extrapolation of solar flare/CME behavior to the more distant stars should be 
appropriate. The Sun has a low overall rate of mass loss (M 0 = 2 x 10~ 14 M 0 yr _1 ) i jof which 


trans ient mass loss events comprise a minor component (typically <10%; 


Howard et al. 


1985 ). However, if mass loss due to CMEs scales with flare occurrence, stars with a high 


flaring rate ought to also s ho w an enhanced rate of mass lost due to transient events 


like CMEs. Recently, 


Aarnio et al. 


(2012) used an empirical scaling between solar flare 


X-ray energy and associated CME mass and extrapolated to pre-main sequence stars with 
measured X-ray flare energies to deduce mass loss rates of putative CMEs associated with 
the stellar flares. They inferred high levels of C ME-related mass loss in these solar-type 
pre-main sequence stars of 10~ 12 -10 -9 M 0 yr -1 . 


Drake et al. 


(2013) also used an empirical 


relationship between solar flare X-ray energy and associated CME mass to extrapolate 
to observed stellar coronal X-ray luminosities and investigate the radiative and kinetic 
energy requirements, assuming that the stellar coronal emissions are produced by flares 
whose occurrence has a power-law dependence on flare energy. They likewise find a large 
CME mass loss rate. The two studies were limited to the X-ray spectral region due to the 
large databases of solar flares measured at X-ray wavelengths, i.e. the GOES 1-8 A band, 
and the fairly large databases of stellar X-ray flares which have accumulated over the few 
decades of sensitive X-ray astronomical telescopes. However, solar an d ste ll ar flare studies 


do take place in other wavelength regions, notably optical (see 


for studies of M dwarf flares and the recent results of 


Hawley & Pettersenl (Il99l| ) 


Maehara et al. 


( 20121 ) for studies of 


energetic flares on solar-like stars). Flare frequency distributions have long been constructed 
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for optical wavelength flare emissions fjLacv et al 


19761 ). and a different method must be 


employed to investigate possible transient mass loss associated with stellar flares studied 
in these wavelength regions. In addition, recent solar flare studies have demonstrated 
that the X-ray band is a minor contributor to th e total radiated fla r e energy integrated 


over all wavelengths, the bolometric flare energy flWoods et al. 


2006; 


Kretzschmar 


2011 ). 


Multiwavelength stellar flare studies have also s’ 


dwarfed by optical flare emissions (Hawley et al 


rown that coronal energy flare losses are 


1995). Taking these factors into account 


requires a better way to intercompare solar and stellar flare energies for a more robust 
determination of the flare bolometric energy. Extrapolating solar flare energy and transient 
mass loss to the stellar regime also needs a way to connect flares and CMEs that goes 
beyond empirical correlations, to avoid the pitfall of “big” scaling with “big” without a 
basis for the relation. In this paper we explore stars with high flaring rates using studies 
from disparate wavelength regions and probe the implications for enhanced stellar mass-loss 
rates using a physically motivated relationship between flares and CMEs. 


2. Flare Energy Partition 


Although there have been many more multi-wavelength studies of solar flares than 
their stellar counterparts, a synthesis of recent results demonstrates similar patterns of 
energy partition for the coronal plasma and optical continuum emission in both solar and 
stellar flares. Table [1] explains a number of energy contributions discussed in the paper. 


The existence of solar white-light flares has been known for years, yet it was only with 


the study of 


Woods et al 


( 20061 ) that a solar flare was detected in an equivalent manner to 


stellar observations, e.g. in integrated solar light (Total Solar Irradiance or TSI). In that 
study the peak intensity of the flare in TSI was only 268 ppm above the non-flare emission, 
compared with a factor of 570 increase in the Geostationary Operational Environment 
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Satellite (GOES) 1-8 A band. The flare in the X-ray band, while exhibiting a much larger 
peak flux enhancement, contribu ted o nly ~l%to the bolometric radiated flare energy: 


Eboi/Ex ~105. The study of 


Kretzschmar (20.11) extended this result to fainter solar flares 


by shifting and co-adding the TSI signal corresponding to the times of GOES flares of 
different peak intensities. They found a range of values for TSI/GOES energy ratios, from 
90±10 for the strongest flares, to 330±130 for the flares that are up to 400 times weaker. 
The coronal flare radiation occurs over a wider wavelength range than the GOES band, 
however, and integrating the radiative contribution of coronal plasma over a wider band to 
encompass most of this reve als that the total s olar coronal flare r adiated energy is 20 % of 


the bolometric flare energy flEmslie et ahl 


2012 ). E C()r /= 0.2. iKretzschmarl ( 2011 ) also 


demonstrated that the white light solar flare emission from these disparate flare classes has 
a continuum shape consistent with that of a black-body spectrum at about 9000 K, which 
contributes ~70% of the total radiated energy. These solar flare studies reveal that the two 
main contributors to the bolometric flare energy are the optical continuum emission and 
soft X-ray radiation. 

Multi-wavelength studies of stellar flares are fewer in number, yet reveal basic 


similarities as to the energy partition. The radiative energy 
in the UV-optical wavelength range has been established by 


Dudget for fl ares on M dwarfs 


Hawley Jz Pettersen 


19 91 ): 


the continuum component has a shape which is characterized by a blackbody with a 
temperature at peak of 9000—10,000 K. This component dominates over line radiation 
when summed over all phases of the flare, carrying 90% of the optical radiated energy; 
-£'cont/-£'o P t = 0.9. Due to the high temperature of this continuum component, the U-band is 
commonly used for photometric monitorin g and is found to contain about 16% of the total 


optical energy flHawlev fe Pettersen 


1991), £u/£ opt =0.16. 


Hawley & Pettersen (119911 ) established the dominance of the blackbody continuum in 

























the optical/ultraviolet part of the stellar flare energy budget. Further observations of the 
optical and X-ray components of a stellar flare measured simultaneously (and with complete 
temporal coverage of the flare in both bands) establish the importance of the coronal flare 
contribution to the bolometric flare radiated energy, assuming that the se two components 


are the major contributors to the flare radiated energy. 


Hawlev et al. 


( 19951 ) studied an 


event on the nearby flare star AD Leo (“EF2”) which had optical and X-ray coverage. The 
energy in the U band filter integrated over all phases of the flare is 2.6 xlO 32 erg (their Table 
2), and using the result above for the fraction of total optical flare energy that appears in 
the U band, gives a total optical flare energy of 1.6xlO 33 erg. The integrated flare energy 
in the EUV bandpass of 65-190 A was 8.2xlO 31 erg (their Table 2). In order to correct this 
for the total coronal flare radiated over a broader bandpass and determine E cor , we use flare 
temperature and emission measure values determined from this flare (peak temperature of 
2x10' K, emission measure of 2.5xlO 51 cm -3 ) and determine the fraction of E cor occurring 
in the 65-190 A(0. 065-0.19 keV) band relative to a wider band of 0.01-10 keV to be 0.12, 
and for the GOES bandpass of 1-8 A relative to this wider energy band, to be 0.29. This 
leads to a total coronal radiated flare energy E cor estimate of 6.8xlO 32 erg. This is to be 
compared against the total optical flare energy of 1.6xlO 33 erg, for a total bolometric flare 
radiated energy E ho i = E cor + E opt of 2.28xlO 33 and a fraction E cor /E bo i=0.3. The two 
dominant components to the bolometric flare energy are the optical/UV emission arising 
from the photosphere/chromosphere and the high energy emission from the corona: this 
implies E opt /E hol =0.7. Then the U band carries a fraction E\j /-E bo i=0.11. Since the energy 
in the hot optical continuum component E cont comprises 90% of the optical energy, this also 
suggests that E cont /E ho i =0.6, similar to the solar flare ratio of 0.7. 


Based on the good agreement between these energy partition values for solar flares and 
the stellar flare discussed above, we assume that these numbers are applicable to all flares 
to estimate the bolometric radiated flare energy for flares occurring on stars of differing 





types. There is a paucity of multi-wavelength flare campaigns on stars of differing types, 


and it is possible that the ratios could vary from flare to flare (or s tar to star), bu 


numbers are a starting point for such investigations. We note that 


Hawley et al. 


these 


(2003) 


found that the optical and ultraviolet emission properties and energy budgets of a sample 
of different types of flares were remarkably similar, which supports this approach. Table [2] 
lists these / = E rad /E bol conversion factors, which estimate the fraction of energy released 
in the band considered relative to the total bolometric radiated flare energy. For stellar 


flares observed with the Kepler mission, an d us ing the tabu 


in the Kepler Instrument Handbook (Van Cleve & Caldwell 


ated Kepler bandpass response 


20091), we compute the fraction 


of the energy from a hot blackbody continuum which appears in this wavelength range, 
using a blackbody temperature of 9000 K and a range of 4000 — 9000 A; the fraction 
f?cont,Kepier/-K C ont=0.26. These energy calculations for solar and stellar cases do not take 
into account the contribution to the total flare energetics of other phenomena known to 
occur along with optical and coronal radiation, particularly the kinetic energy contained 
in accelerated particles, so they do not represent a true flare energy budget, but they are 
used to estimate a bolometric radiated flare energy to place disparate solar and stellar flare 
studies on a more common footing. 


3. Connecting Flares and Coronal Mass Ejections 


3.1. Equipartition between CME Energy and Flare Radiated Energy 


A rough equipartition between flows and ra diation is_a natural consequence if magnetic 


reconnection is the driver of t 
corroborate this expectation. 


re eruptive event ( 


Emslie et ah 


'riest fe Forbes 


2000). Recent solar studies 


(J20Q4|) investigated the detailed energy partition 


in two well-studied flare/CME pairs which had constraints on thermal coronal energy, 
energy in nonthermal electrons and ions, and CME kinetic and potential energy, and 
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Table 1: Different Energy Contributions Discussed in the Paper 


E ra d flare radiated energy in a particular wavelength region 

Eboi flare radiated energy integrated over all wavelengths 

E cont flare radiated energy in the blackbody continuum 

Econt,Kepler flare radiated energy in the blackbody continuum which occurs in the Kepler bandpass 
E opt flare radiated energy in the optical and UV, includes line and continuum radiation 

Eu flare radiated energy in the U filter bandpass 

E cor flare radiated energy from coronal plasma 

Egoes flare radiated energy from coronal plasma in the GOES bandpass 

Ecme total energy in the coronal mass ejection 

Eke,cme kinetic energy in the coronal mass ejection 


energetic interplanetary particles. They concluded that the CME kinetic energy dominated 
over the other energies quantified, by large factors. Based on the later determination 
that total_solar_ radiated flare energies can exceed those measured in the X-ray by large 


factors, 


Emslie et al. 


Emslie et ah 


(120051 ) used refinements to the amount of radiated flare energy in 


(120041 ) to show that although the CME contains the greatest fraction of 


released energy, the uncertainties in the different flare components and increased estimates 
of total radiative energies put the mechanical energy of the CME at a rough equipartition 
with the total radiated energy in the flare. A mor e re cent study of the global energetics of 


a larger sample of solar eruptive events (jEmslie et ah 


2012) confirms this relationship; the 


bolometric radiated flare energy, or the flare energy radiated over all wavelengths from all 
components of the flare, was about one third of the total energy of the mass ejection for the 
38 events studied. The kinetic energy of the CME dominates the mechanical energy, Ecme 

~ E KE , CME - 


For large solar flares, there is a high rate of association between the occurrence of a 
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Table 2: Conversion factors for radiated flare energies in different bandpasses relative to the 
bolometric flare radiated energy 


bandpass 

wavelength range 

value of / = E rad /A b oi 

(Sun) 

value of / = E rad /A bo i 

(active stars) 

GOES 

1-8 A 

0.01 

0.06 

SXR 

0.01-10 keV 

0.2 

0.3 

hot blackbody 1400-10000 A 

0.7 

0.6 

U 

3000-4300 A 


0.11 

Kepler 

4000-9000 A 


0.16+ 


1 assuming T B b =9000 K 


flare and an accompanying coronal mass ejection flYashiro et al. 


2006). Because the solar 


flare model appears to be validated in active stars, we expect that the condition of rough 
equipartition between the energy in the CME and the bolometric flare radiated energy 
in solar eruptive events should also hold in active stars, and use this physical model to 
explore the cumulative effect of transient mass loss in high flaring rate stars. Then, if 
Ecme ~ Eke,cme = 1/2 Mv 2 and e = E hol /E cme , we constrain the mass of the ejection from 
the energy in the flare, given an assumption about the velocity of the ejected material. The 
velocity of stellar CME is unknown, but is taken to be the escape velocity of the star. The 
energy partition established in the previous section is used to correct the observed radiated 
flare energy A rad to its bolometric value by E bo i = E rad //. 


Flare energy distributions are commonly made for both the Sun and other stars; they 
are generally described by a power-law distribution of the form 


dN 

dE vad 


AT E, 


rad 


(i) 


where dN/dE iad gives the differential number of flares occurring per unit time per unit 
energy, and E rad refers to the measurement of flare energies in a particular band of the 
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electromagnetic spectrum. The power-law is generally considered valid over the extent 
-Erad.mm to E rad max . Measured values of a in stars and the Sun range from about 1.7 to 2.3 
(Giidel 2007). 


By expressing the number of CMEs per unit time per unit mass range as dN/dM , the 
total mass lost per unit time due to CMEs is determined using 


MqmE — 


I'M cme 
’ M C MB ,m: 


M, 


CME 


dN 

Jm 


dM 


( 2 ) 


with the limits of integration corresponding to the minimum and maximum CME mass. 
Assuming that the equipartition between CME kinetic energy and total radiated flare 
energy applies to most flares and associated CMEs, 


or Mcme — 
and 4, 


-Eke, cme — —— (3) 

AfcMEV 2 /2 = — l -y- , (4) 

2-E ra d/(e/u 2 ), the total mass lost in CMEs per unit time is, using equations 2, 3 


AfcME — 


^rad.max 2E rad dN 


dE\ 


rad 


-^radp 


(5) 


v 2 fe dE rSLd 

We have assumed a constant linear speed v for all CMEs and we have related the limits 
of integration of mass in the first equation to the corresponding minimum and maximum 
flare energy of the associated flare. Inserting the relationship from equation Q] for the flare 
frequency distribution with energy leads to 


AfcME — 


ad,max 2 J-i] 


rad 


KeE dE ra/ i. 


( 6 ) 


^-Brad.min ^ ^ t 

It is intuitively clear that the cumulative transient mass loss rate should be related to the 


total flare rate. This, it turns out, will allow us to simplify equation [6] The total flare rate 










12 


per unit time can be expressed as 


r *-^'rad 


N tot = 


dN 


-^Tad 


dE, 


dE } 


rad 


rad 


rad ,mc 

IE A 

-'-^r ad, min 

K 


K e E~£ dE rad , 


1 — a 


( rpl—a _ pi —a \ 

\^rad,max ^radjinin/ 


(7) 

( 8 ) 

(9) 


Evaluating the integral in equation 6 from a minimum flare energy T radim i n to maximum 
flare energy T rad>max , and using the result of equation 9 for the normalization, the estimated 
total mass lost due to CMEs per unit time is then 

2 N tot (1 - «) (Pad“max “ ^rad^min) 


Mq.me = 


( 10 ) 


v 2 ef (2 - a) (E 1 -? - E 1 ^ ■ ) 

where N tot is the number of flares per unit time occurring between E rad min and £ raf | mM , 
and E rad is the radiated energy in the particular wavelength band being considered. 

We recast equation [TO] in a format which emphasizes the important variables; setting 
R = E ra d,max/T ra d,min as the ratio of the maximum to minimum flare energies observed, we see 


2 N 

Mqme = — — -y-E iad min G(a, R) 
v 2 ef 


( 11 ) 


where 


G(a,R) = 


1 — a (R 


?2— a. 


( 12 ) 


2-a V-R 1- ”- 1, 

The value of a for a flare frequency distribution determine how much weight is put on 
the largest and smallest flares; large values of a give more weight to T rad , m in while smaller 
values of a do the opposite. Because observed values of a lie in a small range, G varies 
more due to R than with a. Table [3] lists values of N tot , a , T radimax and T rad>min from some 
flare frequency distributions, which together with flare conversion factors in Table [2] can be 
used to calculate flare-associated transient mass loss. 








Table 3. Calculation of Transient Mass Loss Using Published Flare Energy Distributions 


Object 

Vesc 

km s —1 

Atot 

flares star -1 yr —1 

-^rad,min 

erg 

ad, max 

erg 

band 

Q 

A/cme 

(M©/yr) 

A^CME, Aarnio 

(M©/yr) 

ATCME, Drake 

(M©/yr) 

Sun 1 

620 

- Solar Mass Stars - 

200 5.6xl0 27 

5.6xl0 30 

GOES 

1.79±0.05 

4xl0 -16 

8xl0 -16 

6xl0 -16 

Sun 2 

620 

1577 

5.6xl0 27 

l.lxlO 30 

GOES 

2.03±0.09 

lxlO -15 

5xl0 _ls 

4xl0” 15 

young suns in Orion 3 

440 

51 

2xl0 34 

9xl0 36 

SXR 

1.66 

2xl0 -11 

4xl0 -13 

2xl0 -13 

EK Dra 4 

670 

1550 

3xl0 33 

3xl0 34 

SXR 

2.08±0.34 

9xl0 -12 

2xl0 -12 

9xl0” 13 

superflaring sun-like stars 5 

620 

3.5 

5xl0 34 

2xl0 36 

Kepler 

2.3±0.3 

8xl0 -13 

3xl0 -14 

2xl0 -14 

young low mass stars 8 

214 

- Low Mass 

9 

: Stars - 

2xl0 34 

3xl0 35 

SXR 

2.2±0.2 

3xl0 -12 

3xl0 -14 

2xl0 -14 

AD Leo 6 

630 

2557 

6xl0 31 

2xl0 33 

SXR 

2.02±0.28 

5xl0 -13 

3xl0 -13 

2xl0 -13 

AD Leo 7 

630 

3790 

2xl0 30 

3xl0 31 

U 

1.82±0.27 

6xl0 -14 

8xl0 -14 

6xl0 -14 

EV Lac 6 

610 

1461 

8xl0 31 

2xl0 33 

SXR 

1.76±0.33 

4xl0 -13 

2xl0 -13 

lxlO -13 

EV Lac 7 

610 

1790. 

lxlO 31 

5xl0 31 

U 

1.69±0.11 

9xl0 -14 

9xl0 -14 

6xl0 -14 

inactive early M dwarfs 9 

630 

26 

8xl0 29 

10 31 

U 

1.97±0.3 

10 —16 

3xl0 -16 

2xl0 -16 

inactive mid-M dwarfs 9 

630 

50 

4xl0 29 

4xl0 32 

U 

1.20±0.1 

4xl0 -15 

2xl0 -15 

2xl0 -15 


1 Taken from 

2 Taken from 

3 Taken from 

4 Taken from 

5 Taken from 


Yashiroetal 


( 20068 ). 


Veronig et al 


( 20021 ) . 


Wolk et al 


( 20051 ). 


Audard et al, 


( 19991) . 


Maehara et al 


( 20121 ) . 














































6 Taken 
7 Taken 

8 Taken 

9 Taken 


from 

from 

from 

from 


Audard ct al 


i 20001) . 


Lacy et al. 


(1970). 


Caramazza et al 


i 20071') . 


Hilton 20111') . 
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3.2. Empirical Relationship between CME Mass and Flare X-ray Energy 


Both 


Aarnio et al. 


020121 ) and 


Drake et al. 


(120131 ) derived empirical relationships 


between solar CME mass and solar flare GOES energies, with a relationship of the form 


AfcME = KmE , 


7 


M-^GOES > 


(13) 


where 7 is a parameter with an empir ical value of ~0.6 in both studies. Although it 


is not explicitly stated in 


Aarnio et al. 


(120121) . E g in both cases appea rs to be the flare 


energy in the GOES bandpass: the figures in both 


Aarnio et al. 

(2011 

) and 

Aarnio et al. 


(20121) use GOES bandpass flare energies. We use these relationships to estimate a 


cumulative mass loss rate if the same empirical relationship holds to the classes of 
flaring stars considered below, separated from solar flare energies by orders of magnitude. 
Because these stellar flare compilations are not generally made in the GOES bandpass, 
we need to apply the flare conversion factors appropriate for the bandpass in which 
stellar flares are compiled. If we denote this radiated energy bandpass E ia d with a con¬ 
version factor / ra d to relate to the bolometric flare energy, then equation [13] can be re-written 


Mgme = K'E 1 


rad > 


(14) 


where 


K m — Km 


( /goes\ 


(15) 


\ /rad / 

Then the cumulative mass loss rate due to CMEs can be expressed according to equation [2] 
and converting the relationship dN/dM into dN/dE raj & using the relationship between CME 
mass and flare energy, 


AfcME — 


N tot (l-a)K M [ E ^’° 


r»l— ex. 


<£,.d ,max 


pi—ct \ 

■^radjinin / ^^rad.min 


K T dE 


rad 


(16) 
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or, 


A/cme — 


AW 1 - a)K' R 


M 


1—a 


J 7 - a + 1 


/ a+lN 

V^rad.max 


rad, min/ 


(17) 


(E f - E- r . ) 7 

V rad,max rad, min/ 1 

Similar to equation ITT1 with R = E ra d,max/ E radjTnin , we express this so that the dependencies 
are clearer: 

Mqme = Atot K M E? ad m - m G 2 ( , y, o>, R ) (18) 


with 


^2(7, a, -R) — 


1 — 0 ! /R 7 -“ +1 -l 


7 — a + 1 \ R 1 a — 1 


(19) 


Using the parameters lV tot , a, -EVad.max and U ra d,min listed in Table [3l and the conversion 


factors for the GOES band and the appropri ate wave lengt 


lost due to CMEs using the relationships of 


Aarnio et ah 


are also calculated. The values for Km and 7 from the 
K m = (2.7 ± 1 . 2 ) x 10 ~ 3 , 7 = 0.63 ± 0.04, and for 


1 region in Table [2_t he mass 


( 201 l|) and 


Aarnio et a 


Drake et al. 


Drake et ah 


(2013) 


( 2011 ) study are 


( 2013 ) the parameters from 


their empirical fit s are Km = 0.03(0.01 — 0.1), 7 = 0.59 ± 0.02. The empirical equations 


Aarnio et al. and 


Drake et al. 


establish for solar CME masses based on solar GOES flares 


return very similar values of mass for a given flare energy, so their implied Mcme should 


be very similar. However the uncertainty on their fit parameters (particular 
Km) is large, leading to a spread of ~ 150-300 in the mass calculated from 
equation 1 , and a spread of 600-2400 in the mass calculated from Aarnio et al. 


y the consta nt 


Drake et al 


’s equation 2 . 


4. Evaluation and Application 

We have argued that we understand the distribution of flares as a function of energy. 
Combining this with the assumption of a common energy partition between flare energy 
in a given bandpass to the total radiated energy, and that a single constant speed is 
appropriate to be applied to all CMEs, we estimate the amount of mass lost due to CMEs 
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in a variety of flaring stars using a physically motivated rationale for connecting flares 
and CMEs. This method can also be applied to the Sun, where several calculations of 
flare energy distributions have been made, and where the total mass lost in CMEs can 
be determined through observations. The flare energy partitions also allow the use of an 
empirical relationship between solar CME mass and solar flare energy to be extrapolated 
over the several orders of magnitude difference between solar and stellar flares and applied 
to stellar flares observed in a range of wavelength regions. 


We take a constant value of CME speed equal to the escape velocity of the star (for 


th e Sun, this is pa 620 


m 


Aarnio et ah 


km s , so not that different from the average flare-associated speed 


(20 1111 ). We take e = 1 to evaluate the case for an equipartition between 


flare bolometric radi ated e nergy and CME kinetic energy, although we note that the e = 0.3 


determination of 
factor of pa3. 


Emslie et al. 


( 20121 ) would suggest a higher CME mass-loss rate by a 


Table [3] lists parameters for flare frequency distributions from studies of several different 
types of stars, and the estimated total mass lost due to CMEs which might be associated 
with them using the above two formulations (equations [TO] and [171) . 


For the flare energy distributions of stars discussed in the following subsections, we also 
estimated the total mass lost using the same energy ranges for all stellar types considered. 
This facilitates comparison of the technique usi ng different stars and flares characterized 


over different wavelength ranges. 


Yashiro et al. 


(2006) showed that above an X-ray (1-8 


A) flare energy of 5xl0 29 erg, all solar flares in their study had CMEs associated with 
them. Below this level the flare-CME association on the Sun breaks down, for reasons 


whic 

erg 


i are n ot ye t understood. This flare energy level corresponds to an E bol of ~5xl0 


31 


Karpen et al. 


( 20121 ) describe the onset of solar CMEs as being due to the start of fast 


reconnection in a flare current sheet, which triggers explosive energy release and ejection. 
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They comment that the available magnetic free energy determines whether a system 
ends up erupting, which, assuming that the final CME and flare energies traces the inital 
magnetic energy, would agree with the decrease in total magnetic energy available in events 
with decreasing CME and flare energies. We take this lowest bolometric radiated energy at 
which solar flares and CMEs start to exhibit a break in close association as indicative of 
the point where there might not be_enough free energy in the system for the mass to eject 


from the star. As discussed in 


Schriiver et al 


(2012|), the maximum flare energies observed 
on active stars appears to be ~10 37 erg. We use an E bo i imin of 5xl0 31 erg and an Eboi.max of 
10 37 erg, and define i?, bo i = A bo i jmax /.E bol>min and -Rffd = -^boi,max,FFD/-Sboi,min,FFD- Then the 
new lV tot between E bo i imin and E bo i imax is 


Atot = A 


tot,FFD 


E] 


bol.min 


1-a ( ol-ct 


(flw - 1) 


( 20 ) 


, A bo l,min,FFD ) (^FFD 1) 

where A totiF FD is the total number of flares star -1 yr -1 between A b oi, minF FD and A bobmaxFFD 
from the flare frequency distributions (FFDs) listed in Table |3l after converting A ra(Fm in 
and E ra d,max to their bolometric equivalents (A bo i imin FFD = A ra d,mm// and similarly for 
A bo i,maxFFD)- These numbers are listed in Table [U 


4.1. Solar-like Stars 


4-1.1. Sun 


For applica 
distributions of 


ion to the Sun, we take the compilations of solar flare frequency 


Yashiro et al. 


frequency distribution of 


Yashiro et al. 


02006 1 for flares with associated CMEs, and the flare 


Veronig et al. 


02002!) ; parameters are listed in Table |3j In the 


020061) paper, for flare energies below those with a 100% association with 


CMEs, the frequency distribution is corrected by the observed association. Both papers use 
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Tabic 4. Extrapolation of Transient Mass Loss to Common Energy Ranges 1 


Object 

band N to t 

7 ^ flares yr -1 star -1 

A/tot,CME 

M 0 yr -1 

young suns in Orion 

Solar Mass Stars 

SXR 

6 xl0 3 

4xl0 -11 

EK Dra 

SXR 

5x10 s 

4xl0 -11 

superflaring sun-like stars 

i Kepler 

3xl0 5 

10 -11 

young low mass stars 

Low Mass Stars - 

SXR 

5xl0 4 

2 xl0 -11 

AD Leo 

SXR 

10 4 

10 -12 

AD Leo 

U 

2 xl0 3 

10 -12 

EV Lac 

SXR 

6 xl0 3 

7xl0 -12 

EV Lac 

u 

4xl0 3 

10 -11 

inactive early-M dwarfs 

u 

4 

8 xl0 -16 

inactive rnid-M dwarfs 

u 

40 

5xl0 -11 


1 Using parameters (/, a, v esc ) noted in Table [3j and N tot from Table [3] 
adjusted to the common energy range of E bo i min =5xl0 31 erg and 
Ebol.max =10 37 erg. 
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the GOES satellites for full-disk 1-8 A solar X-ray variations, and so we t ake the ratio / to 


be 0 . 01 , based on the total solar irradiance measurements of 


Woods et al. 


these energies to a bolometric radiated flare energy. The solar flare study of 


(2006), to correc 


Veronig et ah 


(120021) compiled distributions for solar flares over a 24 year time interval, without regard 


to CME association. This long time span also covers two solar activity cycles, important 
because of the variation of flare and CME occurrence with phase of the solar cycle. They 
found a power-law of a =2.03±0.09 between 5.6xl0 27 and l.lxlO 30 erg. Their total number 
of flares fitted was 37,851 f lares in 24 years, for 1577 flares yr _1 . These parameters are 


listed in Table [3] 


Yashiro et ah 


( 20061 ) fit a power-law of a =1.79±0.05 between 5.6x10 


27 


and 5.6 xlO 30 erg to solar flares with associated CMEs, with ~10 6 flares/(J m -2 ) above 
their lower energy limit of 2xl0 ~ 3 J rri -2 , or 5.6xlO 27 erg. The time span of their study 
was 10 years, so we estimate N tot in Table [3] to be 200 flares yr _1 . 

These two studies return similar estimates for the cumulative amount of mass lost 


due to CMEs. The two total transient mass loss rates , 4x 


2 % and 6 %, respectively, of the overall M 0 . 


Howard et ah 


-16 


and 10 15 M 0 yr h are 


( 19851 ) showed that the average 


CME mass flux at Earth during a 4 year period (1978-1981) was 5% of the total solar wind 


mass loss. 


Webb (11987 1 demonstrated that the CME mass flux at Earth shows solar cycle 


variations, containing 1 0 % o f the total mass loss rate near solar maximum, and 1 % near 


solar minimum, while 


Jackson & Howard (.1993) demonstrated that 16% of the solar wind 


at the maximum of the Sun ’s ac tivity cycle comes from CME mass loss. The numbers 


derived from the work of 


Veronig et ah 


( 2002 1 ) took into account all flares observed in 


the stated energy range. This overestimates the numbers of CMEs occuring with these 
flares, which suggests that the overestimation should be only a factor of a few, and still 
within the previously observed variation in the total solar mass loss rate occurring due to 
transient mass loss events. This generally good agreement provides additional grounding 
for extrapolating these results to higher stellar flare energies and suggests that if the same 
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processes are occurring in solar and stellar eruptive events, this overprediction of the 
transient mass loss rate might be biased high by only factors of a few. 


4-1-2. The Young Suns of Orion 


Wolketah 


(120051 ) compiled a flare frequency distribution for young solar-like stars in 


the Orion Nebula Cluster observed with the Chandra X-ray Observatory. They computed a 
distribution from the ensemble behavior of selected flares on stars with masses in the range 
0.9-1.2 M 0 at an age of ~ 1MY. Using the spectral parameters for a sample star which had 
temperatures near the median values, we determined conversion factors for the spectral 
energy distribution of the flaring plasma between 0.5-8 keV (the range considered in the 
paper) and 0.01-10 keV, which will encompass the majority of the coronal radiated energy. 
The flare frequency distribution is fitted between -Uo. 5 - 8 keV= 10 34 and 6 xlO 36 ergs, which 
would correspond to a range 2xl0 34 -9xl0 36 erg between 0.01-10 keV. We use / tabulated 
in Table [2] for this energy range to conv ert to the total bolometric flare energy. Masses and 


Wolketah 


(120051 ). and we use the average v esc =440 


radii for these objects are tabulated in 
km s' 1 in our calculations. There were 27 flares above the minimum energy fitted, occurring 
on 25 stars (taken from their Table 6), and the elapsed time of the observation was 660 ks; 
this converts to a flare rate of 51 flares yr -1 star -1 above a minimum flare energy of 10 34 
erg. Table [3] lists a mass-loss rate due to CMEs associated with these flares of 2xl0 -11 M 0 
yr -1 , and Table [4] lists a mass-loss rate due to flares in a common energy range between 
5xl0 31 and 10 37 erg of nearly twice that value. These transient mass loss rates are elevated 
above the current total rate of solar mass loss by factors of ?»1000. 
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4-1.3. EK Dra, a young solar analog 


The nearby (d= 31 pc) GOV star EK Dra has an age estimate of roughly 70 MY 


(Guinan et al 


20031), and is generally considered to be a young solar analog. Because of its 


proximity it has been the subject of much attention . We calculate its es ca pe velocity usi ng 


the stellar mass and radius measurements listed in 


Guinan et al. 


(120031 1. 


Audard et al. 


(1199911 used observations of EK Dra with the Extreme Ultraviolet Explorer (EUVE) to 


explore the rate of coronal flares in EK Dra. The energies determined in that paper 
represent the total coronal radiated energy of the flares from 0.01-10 keV, and so we use 


the / value in Table [2] appropriate for that waveband. 


Audard et al. 


(Il999h fit a power-law 


between flare energies of 3xl0 33 and 3xl0 34 erg, with a flare rate of 5x10 5 flares s 1 
above the minimum energy given above, or a rate of 1550 flares yr _1 . 


The mass loss rate listed in Table [3] for CMEs accompanying these flares is 9xl0 -12 
M 0 yr _1 , fa 450 times the total solar mass loss rate. When considering a common flare 
energy range for mass loss determination, the implied mass loss rate is higher, ~4xl0~ n 
M 0 yr -1 , or a factor of 2000 higher than the present day solar mass loss rate. These rates 
suggest that the Sun might have been able to sustain a much higher rate of mass loss in the 
past. 


4-1.4- Superflaring Sun-like stars 


Recently 


Maehara et al. 


(120121 ) reported on the incidence of flares seen on solar-like 


stars in Kepler data. The bandpass of the Kepler mission is relatively broad (4000-9000 A), 
and we do not know what the spectral energy distribution of these flares is. We assume 
that they are similar to solar and stellar flares in that the optical flares dominate the 
energetics, and that the emission in the Kepler bandpass originates almost entirely from 
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the hot continuum, with the same values of blackbody temperature as seen in solar and 
stellar flares. Table [2] lists the fraction / of the bolometric radiated energy that would 
appear in this bandpass under those assumptions. Further characterization of these stars 
is ongoing, and we do not have precise radii or masses, so we take the solar value of v esc . 
For evaluation of equation [TO] we use the number of flares (183) and the number of stars 
(86) that correspond to the energy range fitted (5xl0 34 -2xl0 36 erg) in the frequency 
distribution (H. Maehara, priv. comm.). The data span 223 days, giving a flare rate of 3.5 
flares yr -1 star -1 . This returns a mass loss rate of 8xl0 -13 for CMEs accompanying the 
flares considered (Table E]), and M of 10 -11 M 0 yr -1 for CMEs associated with a higher 
flare energy range (Table |4|) . 


4.2. Low Mass Stars 


Low mass stars provide a stark contrast with flare behavior from higher mass solar-like 
stars: the factor of ~ three difference in stellar mass and radius is accompanied by a very 
different internal structure. The nearly or completely convective internal structure on low 
mass stars has profound implications for how magnetic fields are generated in the interior, 
and consequent implications for how that magnetic flux emerges and interacts with plasma 
above the stellar surface. As discussed above, detailed studies of individual flares on M 
dwarfs do sh ow agr eement with solar flares. The long timescales for activity decay on M 


dwarfs ( West et al. 


2008) means that if CMEs accompany flares they have the potential to 


act as a significant mass loss process over a large fraction of the star’s life. 
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4-2.1. Young Low Mass Stars in Orion 


Caramazza et al. 


(120071 ) examined the flare frequency distributions of young low-mass 


stars in the Orion nebula cluster from the Chandra X-ray Observatory. They studied 
stars with masses in the range O.1-O.3M 0 , and found no difference in the flare rate for 


by 

Wolk et al. 

(2 

005 

)• 

Caramazza et al. 

(2007) 


(2 0051). We assume that the spectral range over which the energ ies in 


(120071 ) are calculated is 0.5-8 keV (the same as in 


Wolk et ah 


2005) 


as this is not explicitly stated in the paper. We also use the same fraction of radiated 
energy in this bandpass to 0.01-10 keV that we used above for the young solar-mass stars 
in Orion, and the / value in Table [2] for the soft X -ray region t o co rrect this energy to a 


bolometric value. Using stellar parameters from 


Getman et al. 


(120051 ) for this sample we 


calculated a media n escape ve locity for_the low mass stars in the Caramazza et al. sample 


to be 214 km s 


-l 


Caramazza et al. 


(120071) studied 165 sources, and found 151 flares to 


which to fit a distribution. Only flares containing 500 counts or more were contained in 
the distribution fit; this corresponds to an energy of 10 34 erg. Their distribution extends 
to a maximum of ~10 4 counts, giving a maximum energy of 2xl0 35 erg. Converting these 
to the energies expected in the 0.01-10 keV range as in §4.1.2 gives 2xl0 34 - 3xl0 35 erg. 
From their Figure 5, the normalization of their cumulative distribution function at a value 
of 500 counts (their minimum flare energy to which to fit a power-law) is 0.2, for 0.2 x151 
flares to consider. Their exposure time was the same as for the study of 


Wolk et al 


(2005) 


of 660ks. Combining these, we arrive at a flare rate of 9 flares yr -1 star -1 . Table [3] lists 
the calculated mass loss accompanying these flares as 3xl0 -12 M 0 yr -1 using equation ITOl 
Table 0] extrapolates to a mass loss between flares spanning 5xl0 31 and 10 3 ' erg of 2xl0 -11 
M© yr -1 . 
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4-2.2. Nearby Active M Dwarfs AD Leo and EV Lac 


The active M dwarfs AD 

Leo and EV Lac 

lave e 

the past ( 

Hawley & Pettersen 

1991; 

Osten et al. 

2010 ) 


considered to be relatively young (age<300 MY; 


Shkolnik ct al. 


1 nearby (d~5pc) and 


2009) with spectral types 


of M3 and M4 for AD Leo and EV Lac, respectively. Because they are so active they have 
been the subject of many flare studies, and published flar e frequency distributions exist fo r 


both of these objects for both coronal and optical flares flLacv et al. 


1976 


Audard et al. 


2000). This offers an opportunity for a consistency check on the approach taken in the 
current paper. We note that it was AD Leo itself from which the / values were calculated 


above in section 2, and so these numbers should be the most robust for consic 
fl ares. We calculated escape velocities for the two stars using stellar data from 


eration of its 


Reiners et al. 


(1200 91). 


Audard et al. 


(120001) determined the flare frequency distribution for these two stars 


using EUVE observations, and determined coronal flare energy values in the wavelength 
range 0.01-10 keV. We used the value of / listed in Table [2] for this bandpass. They 
determined a flare rate of 7 flares day -1 above a minimum energy of 6xl0 31 erg for AD Leo, 
and 4 flares day -1 above a minimum energy of 8xl0 31 erg for EV Lac. These numbers 
lead to a total number of flares of 2557 and 1461 flares yr~ 3 , and a values of 2.02±0.28, 
1.76±0.33, respectively, for AD Leo and EV Lac. Table [3] lists the values of mass-loss 
calculated for these objects using equation ITUl 5xl0~ 13 M 0 yr -1 for AD Leo and 4xl0 -13 
M 0 yr _1 for EV Lac, and 3xl0 -13 M 0 yr _1 for AD Leo and 2xl0 -13 M 0 yr _1 for EV Lac 
using equation [17l There is good agreement between both M dwarfs as well as both 
techniques. 


For these two M dwarfs, we also made use of the flare compilation of 


Lac y et al. 


dl976|), 


who determined flare frequency distributions for flares observed using U band photometry. 
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We used the value of / for the U band determined above and listed in Table [2] to estimate 
the fraction of the total bolometric radiated flare energy appearing in the U filter ban dpass 


We estimated the minimum and maximum U band flare energies fitted in 


Lacv et al.l ( 19761) 


from the figures for each M dwarf (figures 13 and 14 from that paper); these are listed in 
Table [3] We estimate the total number of flares per year occurring on each star in the 
energy range applicable to their fitted flare frequency distributions (their equation 18) 
by using the a and f3 parameters from their Table 3; for AD Leo this is 3790 flares yr^ 1 
above an E TT min of 2xl0 30 erg and for EV Lac it is 1790 flares yr -1 above an E TT min of 


lxlO 31 erg. Converting between the differential ( a , used here and in 


and cumulative ({3, used in 


Lacy et al. 


Audard et al 


2000 ) 


1976 ) forms of the flare frequency distribution, and 


using the sign conventions appropriate to each, gives a transformation of a = 1 — the 
differential index is listed in Table |3] There is good agreement between the equipartition 
method and empirical method of estimating flare-associated transient mass loss. 


The mass loss rates computed from flare frequency distributions for optical and coronal 
flares listed in Table [3] for each of the two M dwarfs differ from each other by about an 
order of magnitude (5xl0 -13 M 0 yr” 1 , 6xl0” 14 M e yr” 1 for AD Leo and 4xl0~ 13 M 0 
yr” 1 , 9xl0” 14 M 0 yr -1 for EV Lac), ffowever, the range of flare energies considered in the 
coronal flare frequency distribution and the optical flare frequency distributions are quite 
different. When converting to bolometric values, the coronal flare frequency distribution for 
AD Leo covers the range 2xl0 32 —7xl0 33 erg while the U band flare frequency distribution 
covers the range 2xl0 31 — 3xl0 32 erg. For EV Lac the bolometric flare energy ranges are 
3xl0 32 — 7xl0 33 erg (coronal flares) and 9xl0 31 —4xl0 32 erg (U band flares). In both cases, 
the coronal flares are more energetic than the optical flares, and each covers only a small 
range of energy (factor of 20-30 for the case of coronal flares, and factor of 4-15 for optical 
flares), but put together and comparing a values, they indicate a continuation of a similar 
trend over a larger energy region. Combining the two flare frequency distributions into a 
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common energy range allows for a better comparison of the results from the calculations 
done in this paper. This is done in Table [I] where the agreement between the coronal and 
optical flare estimates is much closer, ~1CV 12 M 0 yr -1 for AD Leo for both coronal and 
optical flare frequency distributions, and 7 x 10 12 M 0 yr 4 , 10 11 M 0 yr 1 for EV Lac from 
coronal and optical flare frequency distributions, respectively. 


A different technique attempts to place boundaries on cool stellar mass loss, using the 
detection of astrospheric absorption in the blue-ward wing of Lyman a emission as seen 
at high spectral resolution and using this as an indirect means of detecting cool stellar 


mass loss (Wood et al 


2005). The astrospheric signature is interpreted using heliospheric 


models, part of which implies a mass loss rate for a solar-lik e ste llar wind. T 


EV Lac is one of the targets with an astrospheric detection (IWood et al. 


re flare star 
2005 ). Assuming 


spherical symmetry for a stellar wind and the applicability of the heliospheric model 
applied to the astrospheric detection for this active M dwarf suggests a mass loss rate of 
M = 2xl0~ 14 M 0 yr -1 . The extrapolated rate of transient mass loss for EV Lac, using 
both the optical and coronal flare frequency distributions, is about two orders of magnitude 
higher, according to Table [4j Both mass loss methods rely on extrapolations from the solar 
case, and their inconsistency likely reveals where breakdowns occur between the two. One 
immediate discrepancy is apparent in examining the assumption of spatial distribution: the 
astrospheric models assume a spherically symmetric stellar wind, whereas CMEs originate 
from a particular active region. As the CME expands outward, the motion approximates 
that of a segment of a spherical shell. Thus, in the limit of a high rate of flares/CMES 
encompassing the whole surface, this has the same effect as a steady stellar wind. 


4-2.3. Inactive M dwarfs 


Hilton (2011J) observed flares in the U band on different types of M dwarfs, tabulating 
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the flare frequency distributions on early-, mid-, and late-M dwarf spectral types, for stars 
classified as magnetically active or inactive. A cutoff in characteristic Ha emission level 
was used to discriminate the two activity regimes. Because there are not many studies 
of flare frequency distributions towards magnetically inactive stars, we use the tabulated 
cumulative flare frequency distributions for inactive early- (M0-M2.5) and mid- (M3-M5) 
M dwarfs to estimate the impact of flare-associated transient mass loss for these objects. It 
is important to point out that while these are inactive stars, they are still producing flares 
with similar energies to those of active stars — the main difference is the marked decrease 
in the flare rate. In comparison with the very active M dwarfs described above, these stars 
provide a view of the range of transient mass loss that might be expected from stars with 
lower activity levels. 


The flare frequency distribution is constructed in the same way as for the 


Lacy et al. 


(119761 ) sample, and we use the same treatment as described above for the flare stars AD Leo 


and EV Lac. Flares from multiple stars were used to construct the distributions, and we 
corrected for this effect in determining our N tot . Using the parameters tabulated in Hilton 


(2 0111) ’s Table 4.3, we determine N tot of 26 flares yr -1 star^ 1 above a U band flare energy 
of 8xl0 29 erg for the inactive early M dwarfs and 50 flares yr -1 star^ 1 above a U band flare 
energy of 4xl0 29 for the inactive mid-M dwarfs. Table [3] lists the inferred mass-loss rate if 
coronal mass ejections accompany these flares. Due to the lower flare rate the mass-loss rate 
is correspondingly lower, being 10~ 16 M 0 yr^ 1 and 4xl0 -15 M 0 yr~ x for inactive early- and 
mid-M dwarfs, respectively, using equation [101 an d 3xl0~ 16 M 0 yr -1 for inactive early-M 
dwarfs, 2xl0 -15 M 0 yr -1 for inactive mid-M dwarfs using equation [TTJ These levels are 
similar to those found for solar CME-related mass loss. 


The a value for the inactive mid-M dwarfs is much lower than for flare frequency 
distributions from the other types of stars considered: 1.2 compared to an average of 
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1.94±0.06 for the other distributions listed in Table [3] The CME mass loss rate extrapolated 
to a common energy range tabulated in Table SI gives very different results. This is due in 
large part to the different a values for the flare frequency distributions. An a of 1.2 as for 
the inactive mid-M dwarfs weighs larger flares more heavily, so extrapolating to an energy 
range that encompasses much more energetic flares produces a higher M. Since these stars 
were chosen specifically to have a lower rate of overall magnetic activity as well as flare 
activity, it is not clear if the upper flare energy is applicable. Thus the total mass loss rate 
calculated likely has a higher level of uncertainty associated with it. 


5. Discussion 

5.1. Investigation of Assumptions 

This analysis is predicated on the assumption that the large flares that are observed to 
occur on stars are accompanied by large ejections of mass, and that we can use the observed 
flare frequency distributions to infer the cumulative mass lost from these energetic events. 
Flares manifest themselves in more than one region of the electromagnetic spectrum, 
stemming as they do from structures that span a large vertical extent of the stellar 
atmosphere, from photosphere to corona. We have constructed a physically motivated way 
to connect flares and CMEs, which can be applied to flares observed in different parts of 
the electromagnetic spectrum, by taking into account the fraction of the total bolometric 
energy in that band. The calculated mass loss rates in Table [3] for magnetically active stars 
are large, while in comparison the cumulative mass lost in solar CMEs is a small fraction of 
the total solar mass loss, as expected from observations. 

There are several factors in the method that could affect the hnal calculation: the 
assumption for the flare that E-^ 0 \ = E opt + E cor \ the universality of energy partition from 
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flare to flare; the assumption of equipartition between flare bolometric energy and CME 
kinetic energy (e=l); the error on a, the index of the flare frequency distribution; how 
studies define a flare, which affects lV tot ; the assumption of v — v esc . We discuss those here 
in turn. 


5.1.1. Eboi — £ op t + E, 


Our assumption in tabulating the total bolometric flare energy is that the line and 
continuum emission in the optical/UV, together with the coronal radiation, dominate the 
flare radiated energy budget. The contribution to the radiative flare energy from radio 
and hard X-ray nonthermal radiation will be negligible, due to the radiative inefficiency 


of nonther mal g yrosync 


processes flAschwanden 


irotron and nonthermal bremsstrahlnng compared to collisional 


2002). It is possible that there are other components of flare 


radiation in relatively unexplored 
the radiated flare energy budget. 


wave length regions, 


Christian et al. 


Dut they likely do not dominate 


020031) discussed a flare observed in the 


extreme ultraviolet with a large release of energy in the 60-200 A bandpass, which also 
displayed enhanced continuum emission in the 320-650 A range. The long wavelength 
continuum component contrib uted only at most 10% of the estimated total coronal radiated 


energy. 


Kowalski ct al. 


(120131) described results of a spectral atlas of M dwarf flares which 


showed the commonality of a continuum component appearing at wavelengths longer than 
about 4900 A. While this component becomes energetically more important during the 
gradual phase of the flare, the hot blackbody at shorter wavelengths dominates the optical 
energy. 


Because of the flow of energy through different parts of the atmosphere during a flare, 
it is worth considering whether energy is being double-counted by summing up different 
radiated energy components into a bolometric radiated energy. Radiative hydrodynamic 

























31 


modelling suggests that the black-body which figures prominently in the optical continuum 


emission arises because of nonthermal electrons 


of chromospheric evaporation (Fisher et al 


fitting the lower atmosphere, and theories 


1985) explain the time correspondence between 


signatures which trace accelerated particles (like the impulsive continuum emission in the 
optical, but also radio and nonthermal hard X-ray emission) and signatures of coronal 
plasma heating as the transfer of energy from the accelerated particles to heating plasma to 


coronal temperatures, the Neupert effect flNeupert 


19681) . Backwar ming of coronal X-r ays 


shining on the chromosphere/photosphere is likely to be negligible (Allred et ah 


20061) . 


Additionally, we are not considering the kinetic energy of the accelerated particles in this 
calculation, although other studies of solar and stellar flares have suggested that it could 


dominate over the bolometric flare energies flEmslie et ah 


2012 ; 


Smith et ah 


20051 ). 


5.1.2. Universality of f 


By accounting for the fraction of the total flare radiated energy appearing in a given 
bandpass, flare measurements made in different spectral regions are compared against 
each other and the implication for flare-associated mass loss can be expanded. There are 
relatively fe w ste llar flares which have been observed using a multi wavelength approach. 
Hawlev fe Pettersenl (1991) demonstrated that the optical/UV flare energy budgets of flares 
on M dwarfs appear to be similar to each other, so we are justified in using the same 
fractions for different flaring M dwarfs. The comparison of the flare energy partition values 
between these active stars and the Sun shows generally good agreement, which we take 
as further justification for using the / values derived for flares on different types of stars. 
Further multiwavelength flares observations will be needed to investigate this assumption. 
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5.1.3. Assumption of e—1 


As discussed in §3, there is a reason to suggest a connection between the energetics o 


mass loss and radiation. Solar eruptive events appear to bear this out. 


Ernslie et al. 


( 2012 ) 


described empirical attempts to cha racte rize the energy budget for 38 solar eruptive events, 


and find that EicE,CME/Efl a re,boi ~3. 


Drake et ah 


020131 ) also find a nearly linear relationship 


between the CME kinetic energy and X-ray flare energy for a sample of CME-flare events. 
They restrict their analysis to the GOES energy band, but find that E K e,cme/Egoes ~200 


(their Figure 1), which reduces to EKE,CME/Eboi ~2 after accounting for th e 1% 


of bolometric radiated flare energy appearing in the GOES band. 


Ernslie et al. 


raction 


2005 ) 


noted that considering the uncertainties in the energy estimations, it was possible that the 
energies in the two components are ro ughl y equal. Taking a va lue for e that is 0.5 or 0.3, as 


the solar studies of 


Drake et ah 


(20l3) and 


Ernslie et ah 


(2012), respectively, suggest, would 


lead to mass loss rates that are 2-3 times higher than what is calculated using Equation |TQ1 
H3 or their extrapolated values. 


The association between flares and CMEs must break down at some low energy, where 
there isn’t enough magnetic free energy available for the system to erupt. This has been 
parameterized as the minimum energy at which the solar flare/CME association departs 
from unity, but changing this number downward doesn’t affect the numbers that much. In 
reality, this limit will depend on other factors which will vary from flaring region to region. 
The fact that the observed solar flare/CME relationship departs from unity for smaller 
events demonstrates that there are physical limits on whether reconnection produces an 
eruption which can be observed as a CME. The approach put forward here concentrates on 
the impact of the largest events. 
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5.1.4- Error on a 

The uncertainty on the determination of a leads to a spread in the derived transient 
mass loss rate, but it is a small factor. For most of the published flare frequency distributions 
discussed earlier, varying the index a of the distribution produced only a modest change in 
the derived mass loss rates. These range from as little as a few percent difference to 60% 
difference. This effect is much smaller than the other factors affecting the method. We do 
point out the singularity in equation [TO] for a value of a = 2 , and the singularity in equation 
[17] for a value of a = 1 + 7 (or a ~1.6), both of which span the range of determined a values 
for stellar flares. These equations can be recast to avoid this singularity in the appropriate 
limit (a —> 2 or a —> 1 + 7 ) for large values of R. But this is not necessary because many of 
the studies cited here have only a modest value of R (i.e. the U band compilation of flare 
frequency statistics for EV Lac has a flare energy range R of only 5). 


5.1.5. Definition of a flare 


The various studies considered here used different methods to define a flare, and this 
will impact the results primarily through the rate N tot . The systematic effect these methods 
have on detecting the faint end of the flare distribution will also sway the determination of 
a, but inspection of Table [3] shows a range of a values generally in agreement with other 
studies quantifying the flare frequency distribution. 


Data obtained from optical light curves are binned, by the nature of the data collection 
method; the exposure times set the noise level against which quiescent and flare emission 


is discerned, and thus the integrated energy/fluence of flares. Hilton 1201111 used a series 


of consecutive flux measurements which are several o ab ove a median quiescent flux, for 


light curves with cadences of a few tens of seconds. 


Maehara et ah 


( 2012 ) used a similar 
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method to identify flares, looking at the distribution of brightness changes between pairs 
of consecutive data points, and only looking at the top 1% of this distribution as a way 


to filter out spurious candidates. Since the Kepler light curves used in the 


Maehara et al. 


(2012) study had integration times of n earl y half an hour, this limits the flares to events 


Lacv et ah 


(19761) used a manual “by-eye” inspection of light 


lasting longer than 1 hour, 
curves to identify flares. 

X-ray astronomic al data cons i sts of event lists, which are b i nned 


for flaring variations. 


Wolk et al. 


(120051) and 


Caramazza et al. 


adapatively to search 


(120071 ) used maximum 


likelihood blocks to identify periods of elevated flux (without bias by time binning) then 
required that there be both a rapid rise to the elevated flux level, and_that the elevated 
flux level exceed the quiescent or characteristic level by at least 120%. 


Audard et al. 


(120001 1 


used an adaptive binning of light curves to determine the statistical distribution of counts 
in quiescence; the beginning and ending of flares were identified as times corresponding to 
flux levels separated from the maximum by 2 a. They attempt to correct for overlapping 
flares. 

Solar JSOES light curves were binned in 1 minute cadence for the flare studies. The 


study of 


Veron ig et al. 


(120021) did not correct for background emission in estimating flare 


energies. They also note that their flare definition used four consecutive points with 
increasing flux as part of the determination of the flare rise. 


All of these methods are sensitive to large flares, with the least sensitive only requiring 


a flare duration of an hour for the flare to be detected. 
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5.1.6. Assumption of v = v esc 


In §3.1, it was necessary to assume that the speed of the CME is constant and equal 
to the escape velocity of the star, in order to integrate the cumulative effect of the CMEs 
through the observed flare frequency distribution. Solar CME speeds span a range from 
about_100 km s -1 to 3000 km s -1 , with average flare-associated CME speed of 498 km 


s 1 flAarnio et al. 


20111). CME speeds on the Sun begin at velocities of ~100 km s -1 and 
range up to a few thousand km s -1 . Solar CMEs with speeds less than the solar wind 
speed experience acceleration up to the solar wind speed, and fast speed CMEs decelerate 


to the solar wind speed, which is approximately the escape speed (Webb & Howard 


2012 ). 


For any particular CME we could be off in estimating the mass, by a factor ranging from 
a factor of about 25 lower to 40 higher using the solar CME speed range. But considering 
CMEs in aggregate, the cumulative effect of mass loss should have a value which is closer 
to the ensemble average. If stellar CMEs have larger average speeds, then this biases the 
total mass lost to a lower value. 


Since the empirically derived stellar CME masses using equation [17] and tabulated in 
columns 9 and 10 of Table 0 have a functional dependence on flare energy which is less 
than linear (power-law exponent ~ 0 . 6 ), the weight of large flares is less using these methods 
than for the physically motivated model developed in this paper, where the dependence of 
inferred CME mass on flare energy is linear. This is seen in the M\ arn i 0 and Mq ra ke values 
in Table |3] being smaller than the M values derived above for flare energy distributions 
which go out to large flare energies. Apart from very active stars such as these examples 
(the young suns in Orion, superflaring sun-like stars, and young low mass stars) which 
have E ra d,max >10 35 erg, the agreement between the two approaches is reasonably good. 
Given the many uncertainties involved in making an association between flares and CMEs, 
we consider the general agreement between this method and the empirical method in 
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Aarnio et al.l (2012) and 


Drake et al. 


(2013) to be a good sign. 


5.2. Astrophysical consequences of high mass loss rates 


The mass loss rates suggested by this analysis are l arge for magnetically active flaring 


stars, and have a number of astrophysical consequences. 


Drake et al. 


(120131) considered the 


effect of CMEs associated with stellar flares in a flare-heated corona, and determined a high 
level of kinetic energy required, in excess of the radiative heating input. A nonstandard 
solar model with a high mass-loss rate might be one solution to the faint young Sun paradox 


(Minton & Malh otra 


200 7). The results from Table [3] for the young Suns in the Orion 


Nebula Cluster and the young solar analog EK Dra suggest that flare-associated transient 
mass loss could have been substantial for at least the first ~100 MY of the Sun’s history, 
although the mass loss rate values derived here are lower than what is necessary to change 
the young Sun’s mass significantly. Thus the total mass lost during the active phases of 
these stars is not significant enough to alter their evolution. 

An enhanced mass-loss rate also has implications for the stellar environment. Flares 
which connect to the planet-forming disk could be important in regulating and/or 
influencing planet formation, and coronal mass ejections might also impact the disk in 
young stars, resulting in disk seeding with processed stellar material. This could help 


explain aspects of radionuclide production in protoplanetary disks (IF eigelson et al. 


2002 ) 


The transient mass loss acts as an enhanced stellar wind, and might influence the planetary 


dynamo (IHevner et al. 


2012). CMEs might also play a role in removing material from 


ebris disks and p ro during the recently o 


(lOsten et al. 


2013]). 


Coleman fe Worden 


Dserved dramatic IR variability, as explained by 


19761 ) calculated that for a stellar mass loss of 


about 6 x 10 12 M 0 yr 1 from M dwarfs, a significant amount of ISM material would be 
supplied by M dwarfs. However, even with this high rate of mass-loss, and the abundance of 
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M dwarfs in the galaxy, the cu mulative mass lost is much less than that contributed by the 


rarer low-mass red giant stars (de Boer 


2004). At young ages for planets around solar-like 


stars in the habitable zone, the enhanced transient mass loss might lead to an enhanced 
compression of the planetary magnetosphere and exposure of the exoplanet atmosphere 
to ionizing radiation. Around low-mass sta rs, t his pote ntial effect is magnified due to the 


proximity of the habitable zone to the star flKhodachenko et ah 


2007 ). 


Drake et al. 


The reality of such high inferred mass loss rates has been questioned before. 

020131 ) inferred mass-loss rates as high as 10 - ^ yr -1 for active stars, requiring up 


to 10% of the stellar bolometric luminosity. 


Lim & White (119961 ) used upper limits on 


millimeter-wavelength radio fluxes to argue that winds from active stars could be no more 
than 1-2 orders of magnitude higher that the present-day solar mass-loss rate, in order 
that the observed centimeter-wavelength (and longer) radio flares not be obscured by the 
optically thick wind. Coronal magnetic topologies on both pre-main sequence stars and 


more evo 


ved acti ve main sequence stars show evidence of a complex geometry compared to 


the Sun (IDonati fc Landstreet 


2009[), and a reduction in the fraction of open versus closed 


field lines might have a concomi tant d e crease in the amount of mass loss from a coronal 


wind. On the other hand, 


Cohen et al. 


(120091 ) demonstrated that polar spots on active stars 


could enhance the mass loss rate from a steady wind, albeit in an asymmetric manner. 


There is a complex magnetic topology involved in active regions on stars, and the 
specific configuration may provide an explanation as to whether a particular reconnection 
event could result in an ejection. Strong overlying fields might participate in the 
reconnection, or alternatively provide a barrier to eruption. On stars with stronger surface 
fields and denser coronal plasma than the Sun this could be difficult to achieve. On this 
point we ought to be able to gain insight from the Sun: a case study of AR 12192, one of 
the largest active regions seen in decades, which produced several energetic flares in the 
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October-November 2014 time frame without accompanying CMEs, can reveal the detailed 
pr ocesses for noneruptive ene rge tic fl ares. Preliminary work on failed solar eruptive events 


by 


Chen et ahl (2013) and 


Sun et ah 


(1201511 does suggest that overlying arcades may prevent 


the eruption of a filament. In the end, only direct observations of CMEs in other stars (or 
more stringent constraints) can address whether and how the above topics are influenced 
by transient mass loss. In this regard, observational signatures uniquely associated with the 
CME which are detectable in nearby magnetically active stars provide a path forward. The 
radio signature of coronal mass ejections on the Sun is the so-called type II burst, which has 
a distinct signature of a slowly drifting radio burst. The type II burst arises from plasma 
emission in the tenuous outer atmosphere of the star, produced as the result of MHD shocks 
from the passage of the fast-moving CME through the atmosphere. From a study of the 
dynamic spectrum of such a slowly drifting burst, the velocity of the shock can be extracted 
(providing an upper limit on the CME speed) and with an assumption of equipartition, 
the mass of the CME derived. A direct measurement of a slowly drifting stellar metric 
radio burst would be another step in validating the solar/stcllar connection with respect to 
the flare process, and would establish a relatively unambiguous observational technique to 
study the influence of transient mass loss in active cool stars. 


6. Summary and Conclusions 

We established an energy partition for stellar flares, which allows for an intercomparison 
of flares observed in different regions of the electromagnetic spectrum. Using this energy 
partition, and applying an equipartition between CME kinetic energy and flare bolometric 
energy, we investigated the implication of CMEs associated with these flares for a variety 
of different kinds of stars having measured flare frequency distributions. We relate the 
cumulative transient mass loss to the total flare rate, a point which has not been emphasized 
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in previous studies. While empirical or scaling relations have been used in the past to 
connect flares and mass ejections on stars, the recent details of energy budgets in solar 
eruptive events allows for a physical explanation for connecting the two processes. We 
also used empirical correlations established between solar flare X-ray energy and solar 
CME mass to derive transient mass loss rates. The implied rate of transient mass loss, for 
both methods, is generally large. Our general method also allows for an intercomparison 
of flare frequency distributions for the same object taken in different wavebands; the 
comparison done for two active M dwarfs shows good agreement. There are several 
important consequences to an increase in the transient mass loss from magnetically active 
stars, and observational constraints on the presence of stellar CMEs and their relation to 
stellar flares would be invaluable to evaluating their impact. Just as valuable would be a 
refutation of the solar-stellar connection in this aspect of magnetic reconnection. 

A comparison of implications for flare frequency distributions of the same star at 
two different wavelength regions reveals consistency, which is reassuring. This means that 
the distributions are sampling the same intrinsic flare distribution. The flare frequency 
distribution is commonly used in coronal studies to examine the coronal heating hypothesis 
(suggested for cases where the differential a exceeds 2), but their consistency between 
optical and coronal means that such studies are able to be carried out with optical 
telescopes. The results of the energy partition calculations in Table [2] are used to relate 
the amount of radiated flare energy in one bandpass to the total bolometric amount. In 
particular, the long time baselines allowed by space missions like Kepler, K2, and TESS can 
inform coronal studies, especially in situations of infrequent but large flare energy releases. 
Such missions will increase the yield of flares studied on different classes of stars, at the 
same time expanding the stellar types for which transient mass loss can be investigated. 

RAO acknowledges the useful discussions which took place at the International Space 


Science Institute workshop on Energy Partition in Solar and Stellar Flares for helping 
to frame elements of the study. SJW was supported by NASA contract NAS8-03060 
(Chandra). 
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